Interlayer tunneling in double-layer quantum Hall pseudo-ferromagnets 
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We show that the interlayer tunneling I-V in double-layer quantum Hall states displays a rich 
behavior which depends on the relative magnitude of sample size, voltage length scale, current 
screening, disorder and thermal lengths. For weak tunneling, we predict a negative differential 
conductance of a power-law shape crossing over to a sharp zero-bias peak. An in-plane magnetic field 
splits this zero-bias peak, leading instead to a "derivative" feature at Vb [B\\) — 2nh~vB\\d/ ecj>o, which 
gives a direct measure of the dispersion of the Goldstone mode corresponding to the spontaneous 
symmetry breaking of the double-layer Hall state. 
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Through a series of experimental and theoretical pa- 
pers, it is now well-established that a bilayer of two- 
dimensional (2d) electron gases (2DEGs), when the layers 
are sufficiently close, exhibits an incompressible Quan- 
tum Hall (QH) state at a total filling fraction v = 1, , 
even when interlayer tunneling is negligible. Q. The QH 
state is stabilized by the Coulomb exchange interaction, 
which leads to spontaneous interlayer phase coherence, 
with the layer phase difference </>, the "phason" , as the 
corresponding Goldstone mode. Because the layer index 
may also be viewed as a pseudo-spin-1/2 degree of free- 
dom, this state is also known as the QH FerroMagnet 
(QHFM), and the phasons as (pseudo-)magnons. 

The first direct experimental evidence for this mode 
was obtained recently by Spielman, et al. ||]. For bi- 
layers with small separations d, they observed a giant 
and spectacularly narrow zero bias peak in the inter-layer 
tunneling conductance at low T, a feature that naturally 
emerges from our study of the QHFM. We find rich vari- 
ations of the non-linear tunneling conductance G(V) = 
dl/dV as a function of interlayer voltage V, tunneling 
strength Ao, disorder, temperature, and applied paral- 
lel magnetic field B\\ . Two particularly appealing and 
testable predictions of our theory are: (1) a finite-bias 
feature in G(V) at a voltage V B = Tilo(Q = 2irB«d/ (j>o) 
which tracks the collective mode dispersion relation u>(Q) 
as a function of Bu (</>o = hc/e) and (2) the appearance 
of a power-law negative differential resistance at higher 
bias in cleaner samples. For clarity, below we first sum- 
marize different regimes and results, and then give a brief 
sketch of the calculations. 

We focus here on voltages eV below the QH charge gap 
f2, where the only available bulk excitations are phasons, 
described in the absence of impurities by the Lagrangian 
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where v is the pseudo-magnon "sound" velocity, p s is the 
phase stiffness, M I = ^Jhc/eB is the magnetic length, 



Q is the in-plane "magnetic wavevector" Q — 2irBud/(j)o 
due to applied B|| = B\\y, and to = eV/Ti. The only 
other low-energy modes are edge states, which carry the 
in-phase, in-plane current at low temperatures. 

For weak tunneling, characterized by a "naive" dimen- 
sionless coupling 5q = A./p s <C 1, the tunneling current 
density J{V) can be computed in a controlled perturba- 
tion theory (but see below). There are three interesting 
sub-cases. First, for the pure system, we find 
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where, A = CqAo is the tunneling amplitude downward 
renormalized by quantum fluctuations, with Cq < 1 a 
constant defined below. In above Jo = eA/h2ir£ 2 is 
the tunneling current density scale, and Eg_ = hv/£ is 
the ultra-violet energy cutoff associated with the mag- 
netic length I (presumably 0(£l) for djl of 0(1)), T] = 
k B T/2np s is the familiar Kosterlitz-Thouless (KT) ex- 
ponent pi controlling the asymptotic quasi-long range 
order (qlro) of the double-layer QH state. The exponen- 
tial suppression at high-bias is a kinematic consequence 
of the lack of zero momentum excitations in this energy 
range. 

Note that the above I-V displays negative differential 
resistance, and a peak in J(V) as V — > + . The physics of 
this peak is the appearance of excitons, or electron-hole 
pairs from opposite layers as sharply defined low energy 
excitations in the QHFM, whose spectral function is di- 
rectly probed by the tunneling conductance. 

Second, sufficiently strong impurities modify the trans- 
port (and other aspects of the physics), introducing a 
number of new length scales. We discuss their effects 
primarily in the regime in which tunneling remains per- 
turbative (but see below). If disorder is not too weak 
(weak impurities do not significantly modify the behav- 
ior) , significant quenched-in static deformations in <f> ap- 
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pear beyond a length £, removing the constraint of in- 
plane momentum conservation. Then for eV <C k B T \ 
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where naively a — 77, but preliminary investigations sug- 
gest that at low temperature quantum effects can lead to 
an enhanced a < 1. At large bias, eV ^> hv/t; ^> k B T, 
the current decays exponentially (Id ~ exp(— eV^/hv)). 

Third, in the presence of an additional weak B\\ field 
the zero-bias peak of Eq.|] is split by V B — hvQ/e, 
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Eq. |] is valid for weak impurities, Q£ ^> 1, except in a 
narrow region \eV — TivQ\ < fiv/£, in which the singular 
peak is rounded by disorder and strong-coupling physics. 
Consequently wc find that the location of the tunneling 
conductance feature (which is a peak for Q = and is 
of a "derivative" shape form for Q > 0) maps out the 
wavevector (Q) dispersion of the Goldstone mode asso- 
ciated with the transition to the double-layer QH state 
and with the broken U(l) symmetry. 

The naive application of the above perturbative re- 
sults, I uns = J(V)L 2 (for a square L x L sample), can 
break down for several reasons. First, in the experimen- 
tally relevant geometry with contacts at opposite sample 
ends, there is necessarily a variation of the chemical po- 
tential in the plane of the 2DEG in order to maintain 
the in-plane current. Thus the V in Eqs. |2j-|] becomes 
spatially dependent, and moreover, for sufficiently large 
samples most of the voltage drop is in-plane. Second, and 
more interestingly, the perturbation theory itself breaks 
down at low bias. This occurs once the voltage scale 
Ly = hv/eV becomes longer than a "penetration depth" 
X(V, T) (provided A < L) defined below. Physically, A 
represents a screening length beyond which the in-plane 
out-of-phase currents become too large to support a uni- 
form tunneling current in the absence of another channel 
(e.g. contacts or edge states) to carry away the trans- 
ferred charge. 

In the pure case, we find the non-equilibrium penetra- 
tion depth X(V,T) w A (A /Ly)(/c B T/eV)''/ 2 . A rather 
idealized experimental geometry for which this predic- 
tion may directly apply is for two parallel contacts to the 
same side (but different layer) of the bilayer. Then the 
transport occurs only within a distance A of the contacts, 
within which the tunneling current density is compara- 
ble to that in Eq. || and therefore the current is given 
by Iscr = LX(V,T)J(V). Note that in the "perturba- 
tive" high-bias regime eV/k B T > [<5 (^/A; B T)] 2/(4 ~' 7) 
(equivalently defined by Ly < X(V, T)), X(V, T) is always 
longer than A oqui i(T) = X (X k B T /hv) 11 /^-^ (a length 
that one would have naively expected from an equilib- 
rium calculation), but reduces to it for a lower bias. For 



large tunneling, 5q > I, the penetration depth is instead 
given by a significantly shorter strong coupling expression 
Ao = I\j2-Kp s / A. For the more experimentally relevant 
geometry with contacts on opposite sides of the sample, 
the tunneling current must be carried primarily along the 
sample edges, as in the quantum Hall regime only the 
edge states can carry the symmetric current needed by 
charge conservation. With £ < 00, screening is reduced 
and the corresponding "dirty" penetration depth be- 
comes X d (V) ~ A (Ao/£)(fc B T/el/)"/ 2 - (L V /$X(V,T), 
provided the self-consistent condition £ < A and new 
physics does not intervene at lengths smaller than A^. 

Additional randomness-induced phenomena can also 
occur at scales beyond £, which corresponds to the length 
at which unbound quenched vortices are present. Al- 
though at low temperatures vortex motion is exponen- 
tially suppressed due to the weakness of thermal acti- 
vation and quantum tunneling, because it is neverthe- 
less finite, it leads to temporal fluctuations of the in- 
tcrlayer phase cj> (which we may crudely call "decoher- 
ence" ) for times longer than the exponentially long time 
scale Tvort associated with vortex motion. Therefore, 
provided any unbound vortices are induced, we expect 
an exponentially diverging zero-bias conductance peak 
dI/dV\v=o ~ r vOTt'l^| 2 °f width ^/r VO rt- Ultimately, this 
divergence with reducing temperature is cut off by meso- 
scopic effects in finite samples - phase "diffusion" due 
to energy and current exchange with the leads, tunnel- 
ing contributions from edge modes, etc. - and by effects 
non-perturbative in the tunneling, to which we now turn. 
Finally, for strong enough disorder even the gauge glass 
order may be destroyed at T = 0. The nature of the re- 
sulting state and the transition to it remain unresolved. 

All the above results assume that J(V) (at least in 
a large region near the contact) is perturbative in the 
tunneling A. At sufficiently low V and T, however, this 
perturbation theory breaks down. We first consider the 
clean limit, in which the "naive" dimensionless coupling 
that controls this expansion is S(T) = 5o(£g/kBT) 2 . The 
true expansion parameter is (Lg = 5(kBT/eV)( 4 ~ n ^ 2 , 
as can be seen from standard sine-Gordon RG analy- 
sis. In the strong intcrlayer tunneling regime, defined by 
eV/ksT < <5 2 /( 4-1 '), currents are screened in the bulk on 
the scale A. For contacts on the same sample edge, tun- 
nel currents run only near this edge and the I-V curve 
measures the conductance of the leads. For contacts on 
opposite sample edges, the in-plane current can be car- 
ried only by the edge states, and the model must be ex- 
tended to incorporate these degrees of freedom. However, 
on general grounds (see below) we expect that the power- 
law I-V, displaying negative differential resistance, below 
this strong coupling voltage, generically crosses over to a 
linear I-V, corresponding to a large, narrow conductance 
peak, as observed in experiments. Q 

We now present the highlights of calculations that lead 
to the tunneling results summarized above, defering most 



2 



of the details Rj to a future publication. We focus on the 
tunneling current density J = ^fe^'h ( sm ~ ^ ~ Q x ])> 
where the angular brackets denote an expectation value 
in the non-equilibrium steady state at time t averaged 
over a thermal ensemble of initial conditions in the far 
past. Note that, in the mesoscopic regime in which 
hv/L > eV,kBT, it is crucial to include the degrees of 
freedom of the contacts in the non-equilibrium ensemble. 

Standard time-dependent perturbation theory, to lead- 
ing nontrivial order in Aq gives the current density 
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where C(r, t) = 9{t) ( [ e #(r,*) 5 e -^(o,o)j ^ calculated us- 
ing the unperturbed Hamiltonian with Ao = 0. 
A straightforward calculation gives 

C(r, t) = C Q e - G ^ r ^ (V G Q(^) - c.c.) , (6) 
where C Q = [cosh(n£ e /2k B T)}~ 2 " 1 » e -W^ % an d 
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arc the "thermal" and "quantum" propagators and 
n B (e) = [e e/fcBT - Eqs. f7a||7t] neglect the quan- 

tization of the total charge across the bilayer capacitor, 
and therefore misses Coulomb blockade effects, which 
are small since for large systems, the charging energy 
e 2 /C = h 2 v 2 /p s L 2 is much less than the ballistic "level 
spacing" hv/L. Aside from the presence of Gth, Eqs. |^-^b| 
are identical to the result of classical perturbation theory 
in the equation of motion following from Eq. [I]. 

Gq(t, t) is independent of T and hence reflects the 
quantum dynamics of the system (for A = 0). By in- 
spection of Eq. ^ only when GQ(r,t) contains a non- 
zero imaginary part is C(r,t) non-zero. Therefore the 
tunneling current is non-zero only for vt > r, reflecting 
causality. Since we are interested in low energy scales, 
(eV, hvQ, k B T <C £g) we expect the tunneling current to 
be dominated by the large r,t behavior of C(r,t), with 
p s t > Psf/v 3> h. Then, except in a very narrow range 
(which contributes negligibly to the current) in which 
< vt — r <C hv/p s , Gq <C 1, leading to 
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The effect of thermal fluctuations are encoded in 
Gth ( r i t) . The tunneling current contains two contri- 
butions. Modes at time scales smaller than the ther- 
mal time, tt = h/k^T give T = contributions with 
e -G th 2. in contrast, for modes with vt > r 3> vtt, 
Gth{r,t) 3> 1, with dominant contribution coming from 



thermally activated modes with Tivq <§C k^T, for which 
n B (e) w ksT/e and the upper cutoff on q given by 
ksT/vh. Computing the integral over q in Eq.f7a| under 
such approximation we find 
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Eq. g must be taken together with an implicit causality 
constraint r < vt, (see EqJ|), so that the argument of 
tanh is always real and of order (vr/ty 1 for t 3> tt- The 
well-known KT theory || restricts the applicability of 
Eq. ^| to r) < 1/4, for which the qlro'ed phase is stable 
to thermally-induced vortices. Finally, for bilayers which 
are not too large, an IR cut-off should be introduced to 
restrict r < L. The behavior of C(r,t) for t > L/v, is 
strongly modified by the presence of the contacts even 
at T = 0. This provides an additional source of "de- 
coherence" to smooth the singularities of the low-bias 
tunneling conductance for eV < hv/L. 

We now apply these results for the Green's functions 
to the analysis of the tunneling current in a macroscopic, 
clean system, with L > Ly, where the results are inde- 
pendent of the contact physics. On physical ground we 
expect that for the most parts our predictions will also 
extend into the mesoscopic regime, L < Ly, with the 
mesoscopic modes (not fully accounted for by the </>-only 
model) only contributing a featureless background to the 
tunneling I-V. Combining Eqs. ^ and [9] with EqJ^, we ar- 
rive at a naive perturbative expression for the tunneling 
current density in an infinite sample, 
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which readily leads to Eqs. 

Based on an equilibrium RG analysis, these expressions 
would naively appear to be valid in the thermodynamic 
limit L — > 00, provided eV/k B T > tfW-'O. The true 
regime of validity is, however, more restrictive. In par- 
ticular, Eq. ^ predicts a non-zero steady-state current 
density in the infinite system. This would require that 
either the in-plane (out-of-phase) current or the charge 
density across the bilayer "capacitor" be infinite. The 
apparent violation of charge conservation signals a break- 
down of the perturbative results for L > A. Technically 
it occurs because even large bias V is not sufficient to 
cut off all long scale perturbative divergences. How- 
ever, the remaining divergences are weaker and there- 
fore the length scale X(V,T) over which perturbation 
theory is valid is longer than the equilibrium screening 
length Ae q uii(r) = A (A fc B T /hv)^/^-^ relevant for low 
bias. Physically, A(V) may be estimated by considering 
a steady state, for which the continuity equation reduces 
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to |p s V 2 = J(V), since the in-plane pseudo-spin super- 
current density j s = fp s V0. The naive, uniform J(V) 
can be maintained only over lengths such that the mean 
4> grows by an amount <C 1. Thus X(V) ^ep s /hJ(V), 
which leads to the result quoted in the introduction. To 
support this interpretation, we note that A(V) reduces 
to A equ ii(T) if the naive boundary of perturbation theory, 
eV/ksT w <J 2 A 4 ~'?) ; is approached. For the case of both 
contacts to the same sample edge (x — 0) and L > X(V), 
we have verified perturbatively in A that the non-uniform 
current density profile is J(x,V) = J(V)F(x/X(V)), 
with F(x) — > for x ^ oo and F(x) ~ 1 — 0(\ 2 ). 

The most significant effect of weak disorder is an ef- 
fective random vector potential a(r) (V<f> — ► V0 — a(r)) 
in Eq. [^,[7) This arises from the binding of charge to 
flux in the QH regime, but has short-range correlations 
since it is a physical field. Thus in the zero tunneling 
limit, for weak disorder we arrive at a 2 + 1-dimensional 
XY model with random vector potential (XYRVP). The 
XYRVP is known to have a stable low T phase with qlro 
and no unbound vortices below a critical strength of dis- 
order. H At some threshold strength of disorder, vortices 
will be induced, and the system is effectively in a gauge 
glass regime. The vortices introduce an asymptotically 
exponential decay in r and t 
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C(r,t) -> C(r,i) e - r/ «e-* /T ™ rt . 
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The characteristic length £ diverges as the threshold dis- 
order strength is approached from above, and is finite 
above it, even at T = 0. Vortex localization implies that 
T vort diverges as T — > 0, so r vort ^> £/v. Indeed, there is 
very likely an intermediate regime < f < Tvort over 
which C(r,i) - t-( 1+Q \ with no rigorous theory for a, 
but strong arguments for a < 1. H 

Below the threshold disorder strength, such that £ = 
oo , the XYRVP model may be used to calculate the 
C(r, t). This gives a weak modification of the rj exponent 
by quenched impurities but leaves the results otherwise 
unchanged. Above threshold, the effects of disorder are 
stronger. Eq. [H]may then be used, leading to the pertur- 
bative result in Eq. |^. As for the pure case, such non-zero 
Jd(V) cannot persist in the thermodynamic limit. Simi- 
lar considerations give the result for the dirty screening 
length quoted in the introduction. 

Two different behaviors are expected for weak disorder 
depending upon sample size. If L < A, the low V diver- 
gence of J(V) is cut off by decoherence as phasons scatter 
off the leads. For large samples (L > A), the perturbative 
results are only valid when Ly < A(V), with system in 
a strong coupling, nonperturbative in A regime for lower 
bias. In this limit, we expect a linear response on general 
hydrodynamic grounds. In the presence of tunneling and 
impurities, the only operative conservation laws are for 
total charge and energy. Hence, the low-frequency, long- 
wavelength equations of motion (in addition to that for 
6) must take the form 



where n,; and [Xi are, respectively, the density and chem- 
ical potential in layer i, 7^ = 7(2^ — 1), jjf — 
o~![j dv/ij are the in-plane currents, and repeated indices 
are summed. Note that Eq. [l| combined with the QH 
condition erf? <C erf/, gives a classical description of edge 
transport. The coefficient 7 thereby provides a hydrody- 
namic definition of the tunneling "conductivity" , closely 
related to the true measured inter-layer conductance. For 
zero tunneling, n% — n-i is conserved, and model F (xy- 
ferromagnet) hydrodynamics should apply. | jlO| Hence, 
for small A, we expect that 7(A) is determined by the 
classical crossover between these two limits, which to our 
knowledge is not understood. A scaling approach which 
assumes a direct crossover suggests (for weak disorder) 
7 ~ (fc s r/n)(^/fc B T) 4 /( 4 -")(A/p s ) 2 /( 4 -''). At low T 
(especially for stronger disorder), where the renormal- 
ized A becomes large already in the quantum regime, 
the crossover to the hydrodynamic limit takes place in at 
least two stages, and is more difficult to estimate. The 
general picture which emerges, however, is that tunnel- 
ing currents are screened to sample boundaries. Even 
at low T, this is expected to lead to linear transport as 
V — » 0, and for particular parameters this limit can be 
analytically studied using Eq. 

Although a detailed analysis of the strong coupling 
crossover is difficult, for very weak disorder, the above 
scaling considerations are corroborated by a complemen- 
tary semiclassical analysis, which treats the kinetic and 
tunneling terms on equal footing. As a test of princi- 
pal, we have considered in detail an idealized pure sys- 
tem with parallel side contacts. We find that the cur- 
rent is then carried by a "half soliton" pinned to the 
contacts. Dissipation arises due to a non-local Caldeira- 
Leggett coupling arising from the leads. This gives a 
linear I-V curve at low bias with a conductance that has 
the same scaling form and magnitude as predicted by the 
hydrodynamic approach above. 
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